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Abstract 

Assuming the lightest neutralino solely composes the cosmic dark matter, we examine the con- 
straints of the CDMS-II and XENON100 dark matter direct searches on the parameter space of the 
MSSM Higgs sector. We find that the current CDMS-II/XENON100 limits can exclude some of 
the parameter space which survive the constraints from the dark matter relic density and various 
collider experiments. We also find that in the currently allowed parameter space, the charged Higgs 
boson is hardly accessible at the LHC for an integrated luminosity of 30 fb~ 4 , while the neutral 
non-SM Higgs bosons (H ,A) may be accessible in some allowed region characterized by a large u. 
The future XENON100 (6000 kg-days exposure) will significantly tighten the parameter space in 
case of nonobservation of dark matter, further shrinking the likelihood of discovering the non-SM 
Higgs bosons at the LHC. 

PACS numbers: 14.80.Da,14.80.Ly 
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Introduction: The existence of non-baryonic cold dark matter (DM) has been estab- 
lished by cosmological observations [1]. Weakly interacting massive particles (WIMPs) are 
the natural candidates of DM, among which the lightest neutralino Xi m the minimal su- 
persymmetric standard model (MSSM) has been most extensively studied 

The most convincing detection for the neutralino DM is the underground direct detection 
experiments like CDMS and XENON, which search for neutralino-nucleon (%iV) scattering 
in a low-background circumstance [3, 41. Recently, both CDMS-II and XENON100 reported 
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their search results |3|,|4|], which immediately stimulated some theoretical works [5|,l6(- On the 
theoretical side, great efforts have been paid to improve the accuracy of the prediction for the 
xN scattering rate 3,18]. For example, it has long been known that the hadronic uncertainty, 
especially the strange quark content in a nucleon, can affect the rate by almost one order of 
magnitude, and is therefore impacting significantly the interpretation of the experimental 
searches for DM [9]. This problem was recently better solved by lattice simulation and it is 
found that the strange quark content is much smaller than previously thought, which leads 
to a significant suppression on the uncertainty [101 ]. 

In light of the above experimental and theoretical progress in DM study, we in this work 
re-investigate the xN scattering and use the recent CDMS-II/XENON limits to constrain 
the MSSM parameter space. Note that unlike most recent studies which try to explain the 
two possible DM events reported by CDMS-II, we use the CDMS-II 90% upper limit on 
the spin-independent (SI) xN scattering cross section. In addition to the direct detection 
limits from CDMS-II and XENON100, we also consider the constraints from the DM relic 
density and various collider experiments. We will first perform a scan over the MSSM 
parameter space by considering these constraints. Then we investigate the xN scattering in 
the surviving parameter space to demonstrate the further constraints of CDMS-II/XENON 
on it. Given the extreme importance of Higgs search at the LHC and the strong correlation 
between the xN scattering and the Higgs sector, our study will be focused on the MSSM 
Higgs sector. 

X-ZV scattering in the MSSM: For the sensitivity in current DM direct detection experi- 
ments, it is sufficient to consider only the SI interactions between x\ an d nucleon (denoted 



by f p for proton and f n for neutron 



111 ]) in calculating the scattering rate. In the MSSM, 



these interactions are induced by exchanging squarks or neutral Higgs bosons at tree level 



111] . For moderately light Higgs bosons, the latter contribution is dominant and f p is 
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approximated by 



11] (similarly for f n ) 



fp- J2 — m vfS + l^ha J2 —™ P , 
m a q 27 ^— ' m„ 
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q=c,b,t 



where denotes the fraction of m p (proton mass) from the light quark q while Jt g = 
1 — 2~2 U d s fr i s the heavy quark contribution through gluon exchange, is the coefficient 
of the effective scalar operator given by 11 | 
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with C standing 



then given by 



or the corresponding Yukawa couplings. The ^-nucleus scattering rate is 
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m x o + rriT 

where rriT is the mass of target nucleus and n p (n n ) is the number of proton (neutron) in the 
target nucleus. 

From the above formulas we can infer in which situation the scattering cross section is 
large. Eq.()2]) indicates that this occurs only when Cs xx and/or Cs qq (S stands for a Higgs 
boson) get enhanced. Since the potential enhancement of Cudd by tan is well known, 
we here only analysis the behavior of Cs xx with the variation of SUSY parameters. For a 
bino-likey^ encountered in this work, this coupling is generated through the bino-higgsino 
mixing 2], so a large Cs xx nee ds a large mixing, which means a small fi. To make this 
statement clearer, one may consider the limit of Mi <C M 2 ,/i with M 1; M 2 and fi denoting 
respectively the masses of bino, wino and higgsino. After diagonalizing the neutralino mass 
matrix perturbatively, one can get [5| 

(4) 
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\x cos 2(3. 



So both couplings become large when [i approaches downward to M\. 

In our numerical calculations for the scattering rate, we considered all the contributions 
known so far, including the QCD correction, SUSY-QCD correction as well as the contri- 
bution from high dimensional operators Hj. Note that the SUSY-QCD corrections are not 
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negligible because they may sizably reduce the scattering rate by suppressing Cs qq |8|, H2J . 
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In our calculations we take f% = 0.023, ffj = 0.034, f£ = 0.019, ffj = 0.041 and 
fr^ = fr^ = 0.020. Note the value of fc a we choose is much smaller than that taken in 



most previous studies. This small value comes from the recent lattice simulation [lOJ], and 
it can reduce the scattering rate significantly. 

Numerical results: We make some assumptions to reduce the number of free parameters 
before our scan. First, we note that the first two generation squarks may be heavier than 



about 400 GeV from the Tevatron experiments [13] , and thus their effects on the scattering 
should be unimportant in the presence of light Higgs bosons. So, for the first two generation 
squarks we fix the soft masses and the trilinear parameters to be 1 TeV. We checked our 
conclusion are not affected by such specific choice. Second, since the third generation squarks 
affect the Higgs sector significantly, we let free all the relevant soft parameters. But to 
simplify our analysis, we assume m^, = and Af, = A t , which is well motivated by 
the mSUGRA model with large tan/3. Third, we note that, although the slepton masses 
do not directly affect the xN scattering rate, they can affect the allowed range of tan j3 
via the muon g — 2. In order to avoid a tight constraint on tan /3, we assume a universal 
soft parameter mg and vary it in our scan. Finally, we use the grand unification relation 
3Mi/5ai = M 2 /a2 = M^jctz for the gaugino masses. 

With the above assumptions, the free parameters remained are scanned in the ranges: 
1 < tan/3 < 80, 80 GeV < m A < 300 GeV, 30 GeV < M 1 < 500 GeV, 100 GeV < 
/i, m,j>, rriQ 3 , mfj 3 < 1 TeV and —3 TeV < A t < 3 TeV. In our scan, we consider the following 
constraints like done in [14|: (1) Direct bounds on sparticle and Higgs masses from LEP and 
Tevatron experiments. (2) LEP II search for Higgs bosons, which includes various channels of 
Higgs boson productions. (3) LEP I and LEP II constraints on the productions of neutralinos 
and charginos. (4) Constraints (2c) from precision electroweak observables plus Rb 16] , and 



also the constraints from B-physics observables such as B — > X s 7, B s — > /i~ 



ur, B + ->■ r + u, 



and the mass difference AM^ and AM S . (5) The muon g — 2 constraint 17j (we require 



the MSSM contribution to explain the deviation at 2a level). (6VWe require Xi to account 
for the WMAP measured dark matter relic density at 2a level 1[. The samples surviving 
the above constraints will be input for the calculation of the xN scattering rate. Note that 
most of the above constraints have been encoded in NMSSMTools jl5[. We extended the 
code to the MSSM case, especially we wrote the code for the xN scattering rate to improve 
the scan efficiency. 
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FIG. 1: The scatter plots for the spin-independent elastic cross section of %iV scattering under 
the constraints of dark matter relic density (2a) and various collider experiments. The '+' points 
(red) are excluded by CDMS II and XENON100 (90% C.L.) limits, the 'x' (blue) would be further 
excluded by XENON100 (6000 kg-days) in case of nonobservation, and the 'o' (green) are beyond 
the XENON100 (6000 kg days) sensitivity. 

To show the sensitivity of the \N scattering rate to the value of Jt s , we plot the surviving 
samples on the plane of the scattering rate versus the neutralino mass with the new lattice 
value fx e = 0.02 and with the old value Jt s = 0.38 (corresponding to E^jy — 64 MeV in 
jsj]). One can see that the new lattice value of fa a gives a much lower scattering rate. In our 
following results we fixed fr a = 0.02. 

Our scan samples are 2 x 10 11 random points over the parameter space, and after con- 
sidering the constraints, about 6 x 10 7 samples can survive. We find that for some survived 
samples the xN scattering rate can be as large as 10~ 42 cm 2 , which is far above the current 



CDMS II/XENON limits [3j. Requiring the scattering rate not exceed the current CDMS 
II/XENON bounds, we find that about 33% (61%) of the survived samples are ruled out for 
7B4 < 300 GeV (200 GeV). Further, if the future XENON100 with 6000 kg-days exposure 
18[ gives null dark matter results, then 90.5% (99.5%) of the survived samples will be ruled 
out for m,A < 300 GeV (200 GeV). So the dark matter direct detection experiments are 
highly complementary to collider experiments in testing the MSSM. 

In Fig. [2] the surviving samples are projected for tan (3 and /x versus the charge Higgs 
boson mass m#+. We see that the regions excluded by the CDMS-II/XENON limits are 
characterized by large tan (5 and small /i. The future XENON100 (6000 kg-days) can further 
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FIG. 2: Same as Fig. [TJ but showing tan /3 and [i versus m H + for fr a = 0.02. 
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shrink the currently allowed regions, and in particular set a bound m H + > 165 GeV. Around 
this lower bound, the value of /i is quite large (~ 1 TeV) and so the MSSM has a fine-tuning 
due to the relation m 2 H+ > m\ and m\ = m\ u + m 2 hd + 2/x 2 

In our following discussions, we only focus on the samples that satisfy the current CDMS- 
II/XENON limits. In Fig. [3] we project the surviving samples in the planes of tan/3-// 
and tan f3 e Q-m H + . Here tan/3 eff = tan/3/(l + A b ) with A b denoting the SUSY radiative 
corrections to bottom quark mass[l2|. As expected, large tan/3 must be accompanied by 
large fi to suppress the the scattering rate, and this tendency becomes more apparent for 
the samples further satisfying the future XENON100 limit. We note that in this case, i.e. 
large tan/3 along with large /z, A b should be large [121 so that tan/3 cff is significantly smaller 
than tan 0. This speculation is verified by Fig. |3]and also by our results for which show 
A;, larger than 30% for tan/3 > 40. 

Implication for LHC Higgs searches: Above results showed that the CDMS-II/XENON 
limits have set upper bounds on tan/3 e g. Since the LHC search for non-SM Higgs boson 
usually needs a large tan/3 e ff to enhance the signal rate 19l-l21|. such upper bounds on 
tan /3 e ff may have important implication on LHC search for non-SM Higgs bosons. 

We first consider the LHC search for the charged Higgs boson, which, for the charged 
Higgs heavier than top quark, mainly utilizes the channel gg/gb — > t[b]H + with H + sub- 
sequently decaying to r + z/ T 19[. In Fig. H]we show the rate of this channel in the allowed 
parameter space, where the model-independent 5a discovery sensitivity is obtained by the 
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FIG. 3: Same as Fig. [21 but showing tan/3 versus and tan/3 e ff versus m H +. 
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ATLAS collaboration for 30 fb _1 integrated luminosity 19] . In the calculation of the signal 
rate, we used the effective lagrangian method to incorporate the important SUSY correc- 
tions. Our results show that for more than 99% of the survived samples, the rate is smaller 
than the discovery sensitivity, which means that the LHC is unlikely to discover H + . Our 
results also indicate that the future XENON100 limits (in case of nonobservation of DM) 
will further tighten the parameter space, making the discovery of H + unlikely even with 
higher luminosity. Note that for the charged Higgs lighter than top quark, the LHC search 
can instead utilize top pair production with one top decay into charged Higgs 19[. Like the 
case of the heavy charged Higgs boson, our results indicate that for more than 99% of the 
survived samples, the signal is below the 5a discovery sensitivity obtained by the ATLAS 
collaboration due to Br(£ — > H + b) < 10~ 2 . The small branching ratio of t — > H + b arises 
from the fact that tan /3 is around 10 for < 150 GeV (see Fig. [TJ) and for such a value 
of tan (3 there is a strong cancellation between different terms in the amplitude of this decay. 
Now we turn to the LHC search for the non-SM neutral Higgs boson H and A, for 



which both the ATLAS and the CMS col 
bbH(A) with H(A) decaying to r leptons 



aborations utilize the channels gg — > H(A) or 
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21] . Unlike the charged Higgs boson search, 



for which a model-independent discovery sensitivity can be obtained, the analysis for the 
search of the neutral Higgs bosons is performed in certain SUSY scenarios. Here we consider 
the m™ ax scenario with the following fixed parameters: Msusy = 1 TeV, M 2 = 200 GeV, 
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FIG. 4: Same as Fig. [3l but showing the LHC search sensitivity for the charged Higgs boson [19j ]. 
Here the LHC sensitivity curve obtained by the ATLAS collaboration corresponds to the 5a 
discovery level, while the exclusion limit from the dark matter direct detection experiments is at 
90% C.L. 



m-a 



800 GeV, and X t = A t — /i cot (3 = 2 TeV. In order to show the /i dependence of the 



constraints, we choose several representative values of \i and scan the rest free parameters 
in the ranges: 1 < tan/3 < 80, 80 GeV < m A , mi < 0.8 TeV and 30 GeV < M x < 500 GeV. 
In Fig. [5] we show the surviving samples on uia versus tan /3 plane together with the LHC 
discovery sensitivity for 30 fb" 1 integrated luminosity. This sensitivity is obtained by the 
CMS collaboration with H/A — ► t + t~ — > p, + jets topology (semi-leptonic final states) |20|, 
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FIG. 5: Same as Fig.[3j but showing the LHC search sensitivity for the non-SM neutral Higgs bosons 



in the m™ ax scenario [20]. Here the LHC sensitivity curve obtained by the CMS collaboration [2C] 
corresponds to the 5<t discovery level, while the exclusion limit from the dark matter direct detection 
experiments is at 90% C.L. 

which is better than that obtained by the ATLAS collaboration with H/A — > t + t~ — > 2£+4u 



topology (full leptonic final states) 19l. l21|. Note that in Fig. [5] the LHC sensitivity curve 
corresponds to the 5cr discovery level, while the exclusion limit from the dark matter direct 
detection experiments is at 90% C.L.. 

A few comments are due regarding the results displayed in Fig. [5j 
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In getting these results we used the package NMSSMTools (version 2.3.1) 15| which 
uses micrOMEGAs (version 2.2) [22| for the calculation of the dark matter relic den- 
sity. But in our calculations we extended the package by including more experimental 
constraints, such as the LEP search for the Higgs bosons and B + — > t + u t , so our 
combined constraint on the parameter space is more stringent. 

The CDMS-II/XENON constraints are sensitive to the value of fi, i.e., as /i gets larger, 
the constraints become weak. The reason for this behavior is that a larger fi will result 
in a smaller Higgsino component in Xi an d hence suppress the Higgs-Xi-X? coupling, 
which will weaken the CDMS-II/XENON constraints. 



The LHC sensitivity for fi = 200 GeV is taken directly from [20J, and for other values 
of fi the curves are obtained by scaling the value of tan (3 so that the production rate of 
the Higgs bosons is same as that for fi = 200 GeV. In doing this we used the package 



FeynHiggs2.7.1 23) to calculate the production rate. Note that the /i parameter affects 



the production rate mainly by changing the H(A)bb coupling through loop correction 



(A;,) which is proportional to fitan(3/MsusY [12]. So the shift of the LHC sensitivity 



curve due to the variation of the fi value is not negligible, as shown in Fig. [51 

Figj5] shows that for fi = 200 GeV no surviving samples can reach the observable 
level, while for larger /i values a small fraction of surviving samples can lie within 
the observable region. Numerically we checked that for fi = 400GeV, 700GeV, ITeV, 
about 8%, 11% and 7% of the surviving samples lie within the observable region, 
respectively (for fi = 1 TeV about one third of these detectable samples can even 
survive the future XENON100 limit). The reason for this behavior is that for a large 
fi, although the LHC sensitivity curve is shifted upward, due to the much weakened 
CDMS-II/XENON constraints, some surviving samples can have quite large tan/3 
values (as shown in Figj3] and Figj5]) so that they can reach the LHC sensitivity. 

About the lower bound of Ma as a function of tan/3, since both the xN scattering 
cross section and the production rate of the Higgs boson are proportional to tan 2 (3 for 
large tan [3, one would naively expect the lower bound curve runs in parallel with the 
LHC sensitivity curve. As shown in Figj5j this is not the case because we considered 
many experimental constraints and not all of them scale as tan 2 (3. 
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Conclusion: We have seen that if the MSSM is the true story, the current limits from dark 
matter and collider experiments already strongly constrain the parameter space, which has 
important implication for the LHC searches for the non-SM SUSY Higgs bosons. It turns out 
that in the currently allowed parameter space, the charged Higgs boson is hardly accessible 
at the LHC for an integrated luminosity of 30 fb _1 , while the neutral non-SM Higgs bosons 
(H,A) may be accessible in some allowed region characterized by a large \i. The future 
XENON100 (6000 kg-days exposure) will significantly tighten the parameter space in case 
of nonobservation of dark matter, further shrinking the likelihood of discovering the non-SM 
Higgs bosons at the LHC. So the interplay of the dark matter direct detection experiments 
and the LHC Higgs searches will allow for a good test of the MSSM. 

Finally we stress that we obtain the above conclusion by choosing a small fas (= 0.02). 
If we choose a large /y s , the scattering rate will be larger so that the limits from the current 
CDMS II/XENON will become more stringent. For example, for fx s = 0.38 taken in 
previous studies [9(, we find that the current CDMS II/XENON constraints are comparable 
with the future XENON100 (6000 kg-days) constraints. We also checked that if we relax 
some assumptions in our scan, e.g., the grand unification relation for the gaugino masses, 
our findings about the LHC Higgs searches remain unchanged. Further, we noticed the 



controversy on the XENON100 detection efficiency [24J. Although our current bounds are 
from CDMS II plus XENON100, the CDMS II plays the dominant role. If we do not include 
the current XENON100 limits, our results almost remain unchanged. 

Note added: After finishing our paper, we notice that the ATLAS collaboration publishes 



an analysis 



or the search sensitivity of the neutral non-SM Higgs bosons via the semi-leptonic 



final states 25] , in which the obtained discovery sensitivity seems better than the CMS result 
shown in Figj5j According to this new ATLAS result, more surviving samples in Fig J5] will 
reach the observable region. So our conclusion about the observability of the neutral Higgs 
bosons will remain unchanged. 
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